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ABSTRACT 

Electrical contacts and metallizations to GaAs solar cells, 
intended for satellite array and other space-based a pp ^ at ^°" ' 
must survive at high temperatures for several minutes und 
soScific mission scenarios. The determination of which 

SKSmzatxons or alloy system that are able < tc > «^st and extreme 
thermal excursions with minimal degradation to solar cel 
nerforlance can be predicted by properly calculated temperature 
constitution phase diagrams. Alternately, the thermodynamic 
behavior of metal-GaAs systems could only be accomplished y 
performing hundreds of controlled reactions and experimental y 
determining compound formation, diffusion couples, eutectic points 
and finally, electrical/mechanical performance degradation. A 
Method Tor calculating a ternary diagram and its three constituent 
binary phase diagrams is briefly outlined and ternary phase 
diagrams^ for three Ga-As-X alloy systems are presented. Free 
energy functions of the liguid and solid phases are approximate y 
the regular solution theory. The three binary systems comprising 
the boundaries of a ternary phase diagram are utilized to calculate 
the binary regular solution parameters. The free energy functions 
for the ternary system liquid and solid phases are then written a 
a cognation of £he binary regular solution pari meters Liquldus 

and solidus boundaries, as functions of temperature and 
composition, for the ternary diagrams are calculated from these 
free energy functions. Phase diagrams calculated using this method 
are presented for the Ga-As-Ge and Ga-As-Ag systems. 

INTRODUCTION 

This research effort was initiated to determine an alloy system 
which is chemically stable when exposed to extreme te^P^ture 
excursions, exceeding 873 K, while functioning as electrical 
contact metallization to GaAs. The theoretical portior i tb 

investigation has been the calculation and utilization or 
temperature constitution phase diagrams to determine suitable : al oy 
compositions, which will be in equilibrium with the (p) or (n)GaAs 
solar cell emitter when the device is thermally stressed. Large 
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experimentally and are readily available i^the opl^l iterator™ 1 [1^ 
blen expTrTm'entally The a h l0y SySte ” S h * Ve 

^as^^nSS 3 ^ 1 s “a h£rT\^ 

Mathematical representations for the liquidus and c n ^H y stams - 
boundaries were derived through an application of fho ® Jl ldus P hase 

s-ssjssirs 1 . 

parameter is a linear function of the absolute temperature 

solubility for a^of theb^Temlnt^^on^u^tr Snd ( th° 
Neumann-Kopp rule is obeyed [ 8 ] . f the 


BINARY ALLOY SYSTEMS 

SSS" is Sr = 

describes the equilibrium between a binary liquid ^ale and 1 ^ 

terminal solid solubility phase Th* phase and a 

equilibrium between a binary liquid ^ohale 0336 • t + reats the 

binary^regular b ^so'^ti^>n^ 0 parameters? alC ^nb^)^^^ 1 ^Q U ^^^F^ a 
were used to generate the boundaJi^ fo^the Ugu^rand^nTbr 3 

obtained ^biues^o^tlr^pinise ^c^ii l^br^unT boundaries X ^ er ^ men b a 


EQUILIBRIUM BETWEEN A BINARY LIQUID 
SOLID SOLUBILITY PHASE 


PHASE AND A 


BINARY 


TERMINAL 


In developing the mathematical 
thermodynamic equilibrium between 
terminal solid solubility phase 
principles were utilized: 


relations to express the 
a binary liquid phase and a 
the following procedure and 


1. 


2 . 


Apply the principle of heterogeneous phase equilibrium. 

Represent the partial molal free energy terms of each of fh G 
elemental components as functions of the standard state free 
energies of the pure elements, the activity coe«Ic!entI, £hj 
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mole fractions of the components, and the absolute 
temperature . 

3 The difference between the standard state free energies of 

the liquid and solid states of each of the elements is 
expressed in terms of their heats of fusion and melting 
temperatures. The difference between the heat capacities of 
the liquid and solid phases is assumed to be zero. 

4. Use the regular solution model of alloying behavior to obtain 
expressions for the activity coefficients of the elemental 
components . 

A straightforward application of these steps results in the 
expressions for the boundaries of the phase diagram. 

The analytical representation for the thermodynamic equilibrium 
between a binary liquid phase and a terminal solid solubility phase 
is given by the two simultaneous equations: 

a//F 1 x(1 - T/T ± ) = RxTxln (N 1S /N 1L ) + S12xN% s - L12 xN?> l (la) , 


and 


lHF 7 x ( 1 - T/T 2 ) = RxTxln(N 2S /N 2L ) + S12xNf s L12xNf L 


(lb) , 


where HF, and HF 2 are the heats of fusion for element 1 and element 
2 respectively; L12 is the liquid phase regular solution parameter 
for alloy 1-2; S12 is the solid phase regular solution parameter 

for alloy 1-2; N lL , N 2L , N ls , and N 2S are th ® m ° 1 \ fr ^ 1 ° R nS .f 
1 and element 2 in the liquid and solid Phases, R is the gas 
constant; T is the system temperature; and T, and T 2 are the melting 
temperatures of the two component elements. 

The regular solution parameters, L12 and S12, are assumed to be 
linear functions of the absolute temperature as shown below: 


L12 = a 


(5T 


(2a) 


and 


S12 = x + 6T 


(2a) , 


where the slope 
of the data to a 


and intercept constants are determined from a fit 
linear function of the absolute temperature. 
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A BINARY COMPOUND 


EQUILIBRIUM BETWEEN A BINARY LIQUID PHASE AND 
PHASE 


The relationship developed for describing composition of the liauid 
boundary which is m equilibrium with an intermediate binary 
compound was developed by following the following steps: Y 

i. The molal free energy of the compound is expressed in terms 
of the compound composition and the partial free energies of 
its solid elemental components. 


2 . 


3 . 


4 . 


5. 


e erogeneous phase equilibrium principles are used to 
describe the partial molal free energies in terms of the 
elemental components of the liquid phase which are in 
equilibrium with the compound. 

The molal free energy of the compound is expressed in terms 
of the absolute temperature, standard state free energies 
e fractions, and the activity coefficients of the ' 

components in the liquid phase. 

The molal free energy of the liquid phase, which has the same 
composition as that of the compound, is described using the 
same variables as in Step 3 above. y 

The heat capacity differences between the compound and the 
liquid are assumed to be zero. The molal free energy 
differences of the compound and the liquid phase are 
expressed in terms of the heat of fusion and temperature of 
the compound. The free energy difference between the liquid 
a ^ 0 SOlld phases at the equilibrium melting temperature is 

z Biro . 


6 . 


Regular solution theory i 
the activity coefficients 
phase . 


s used to define the relations for 
cf each component of the liquid 


The resulting equation which describes 
binary liquid phase and the binary 
relation [9], 


the equilibrium between the 
compound is given by the 


aH c x( 1 T/T c ) - RxTx (N lc xln (N 1C ) + N 2C xln (N zc ) ) 

- RxTx (N lc xln (N 1L ) + N 2C xln(N 2L )) 

+ N lc xL12xN% c + N 2C xL12xNl c < 3 ) 

- N 1c xL12xN% l - N 2C xL12xlf 1L . 


RESULTS FOR BINARY ALLOY SYSTEMS 

Gallium-Germanium 

To demonstrate the validity of this computational method of 
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calculating equilibrium phase diagrams, the gallium-germanium (Ga 
svstem is discussed in detail as an example and compared to 
SiiHcal data from the literature. Equations (la) and (lb) were 

“S. the binary reguiar solution parameters 

Flemental values for the melting temperatures and the heats of 
“Son terl takin from the compilation of Hultgren [10] and are 
listed in Table 1. Temperature and composition values for 
liauidus and solidus boundaries, required for these calculations, 
Irlfrom the experimental efforts of Keck and Broder [11], Greiner 
and Breidt [12], and Thurmond and Kowalchik [HJ. 

Table 1 


Values for the Heats of Fusion and Melting Points of Selected 

Elements __ == _ 


Element 

Heat of Fusion 
(Cal/gm-at) 

Melting Point 

(K) 

Silver (Ag) 
Arsenic (As) 
Gallium (Ga) 
Germanium (Ge) 

2700 

2600 

1335 

8100 

1234 

1090 

303 

1210 


At several temperatures the values for the regular solution 
parameters were calculated. These data points were fit to a linear 
function of the temperature by the method of lea: 3 t squ, ares. 
Calculated values for the slope and intercept for the Ge-Ga . and 
selected binary systems are shown in Table 2. As a cnecK 
the^omputat ions , phase boundaries for the Ge-Ga binary system were 
calculated and compared to experimental values from Hansen [14] • 
Figure 1 shows the liquidus and solidus boundaries ca lc ulated for 
the Ga-Ge system. Ge solubility m Ga is very small and the region 
of terminal solubility is almost coincident with pure Ga h 

comoosition calculated for the eutectic of Ga-Ge is located close 
to P the melting point of pure Ga, which is in agreement with 
experimental results. The calculated and experimental liquidus 
boundary difference is less than 1 atomic Per^nt over th e complete 

calculations ^V'Th^he “"s^me retrograde solubility 

characteristic that has been experimentally observed. Maximum 
SoISbiUW in solid Ge occurs near 650"c for both calculated and 

experimental boundaries. 


GaAs Liquidus Boundary 

Exoerimental values are available for only a few compounds of 
interest [8 9], therefore two sets of computations were completed. 
In^the first set , values for the heat of fusion, melting point and 
composition of GaAs, and the composition and temperature of the 
liquid phase that is in equilibrium with the ^ compound were ^ used to 
calculate the regular solution parameter for the liquid phase. 
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Experimental values for GaAs were taken from References [14 15 191 
The regular solution parameter calculated from these values usinq 
Equation 3 is shown in Table 2. y 


Table 2 


Binary Regular Solution Parameters 



Liquid 

Liquid 

Solid 

Solid 

£3 

System 

Intercept 

Slope 

Intercept 

Slope 

Ref. 

Ag-As 

Ag-Ga 

As-Ga* 

As-Ga** 

As-Ge 

Ga — Ge 

* Calcu 1 at*p 

-9256.9 
-28 ,269 
10,455 
19,264 
-6200.4 
210.23 

H wif-h 

. 12305 
18 . 887 
-11.608 
-11.233 
4 . 688 
-.54553 

134.12 

-25,839 

16,816 

5945.4 

-4 . 1037 
15.932 

-6.0467 

-2.0719 

[16] 

[17] 
[14,15] 
[14,15] 

[18] 
[H-13] 


* Calculated with the heat of fusion = 10,758 cal/gm-at. 
♦♦Calculated with the heat of fusion = 3,613 cal/gm-at. 


The second set of computations involved using the same values for 
the melting point and composition of GaAs and the same liquidus 
boundary to determine values for the heat of fusion of GaAs and the 
regular solution parameter of the liquid phase. The GaAs heat of 
fusion calculation, in addition to those for other binary 

a In e TaiS1. in ^ ^ ^ SlOPe »" d ^ -lue£ 


Table 3 

Heat of Fusion of Selected Binary Compounds 


Compound 

H f Calculated 

H f Experimental 

AsGa 

3613 

10, 578 

As 2 Ge 

7789 


AsGe 

8087 

— 


The heat of fusion value difference between the calculated and 
experimental results is considerable. Similarly, the values 

disoa£itv ^ Tn re 9 ular solution parameter shows a significant 
disparity. in order to determine accuracy of subsequent 
calculations using these values, the liquidus boundaries iere 
computed using the values obtained in the two sets of calculations 

the twZ U ^2 9 n5 a " A i P ? a J e diagram sh °wing the comparison betweeA 
the two sets of calculations are shown in Figure 2. Statistical 

anaiys! 3 of the results shows that the maximum deviation letveen 
the computed and experimental boundaries is less than five atomic 
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percent arsenic. While the thermodynamic values are significantly 
different between the two approaches, the resulting liquidus 
boundaries for the calculated and experimental values are in 
satisfactory agreement. 


TERNARY ALLOY SYSTEMS 

Expressions describing the equilibrium in ternary systems were 
developed for two situations. First, the equilibrium between a 
ternary liquid phase and a ternary terminal solid solubility phase 
is represented by three simultaneous equations. Secondly, the 
expression for the equilibrium between a ternary liquid phase and 
a binary compound was formulated. The thermodynamic behavior is 
approximated by the binary regular solution parameters in both 
instances . 


EQUILIBRIUM BETWEEN A TERNARY LIQUID AND A TERNARY SOLID SOLUBILITY 
PHASE 

The assumptions and procedures for developing the relations 
describing the binary equivalent to this case were applied and the 
following three simultaneous equations were generated. 


ah, x(l - T/T.) = RxTxln(N ls /N 1L ) + L12xN 2[ x( 1 - N 1L ) + L13xN iL x( 1 N 1L ) L23xN 2L xN 1L 

- S12xN xs x[ 1 - N 2S ) - S23xN is x{ 1 - N 2S ) + S13xN ls x N xs , 


(4a) 


aH (2 x ( 1 - T/T 2 ) 


RxTxln {N 2S ) / N 2t H 
S12xN ls x ( 1 - N 2S ) 


L12xN 1l x (1 - N 2L ) 
- S23N Js x( 1 - N 2S ) 


L23xN 3L x( 1 - N 2L ) 
+ Sl3xN ls xN 3St 


- L13xN 1L xN :iL 

(4b) 


and 

aH,,x ( 1 - T/T,) = RxTxln(N iS /N iL ) * L13xN xl x[1 - N XL ) * L23xN 2[ x (1 - N 1L ) 
- S13xN ls x ( 1 - N is ) - S23xN 2S x( 1 - JV 3S ) > S12xN ls xN 2S . 


LI 2xN lL xN 2L 

(40 


Equations (4a) -(4c) were solved to yield a solution, as a function 
of temperature, for the ternary Ga— As— Ge liquidus boundary which is 
in equilibrium with the Ge solid solubility phase. The binary 
regular solution parameters, which are used to estimate the 
alloying behavior of the ternary alloys, are those listed in Table 
2. A 1000K isothermal section of the resulting ternary diagram is 
shown in Figure 3a and includes several tie-lines which connect the 
compositions on the liquidus boundary which are in equilibrium with 
the composition on the Ge solidus boundary. These calculations 
were repeated for the Ag-Ga-As system and an 800K isothermal plot 
of the liquidus and solidus boundaries is shown in Figure 3b. 
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EQUILIBRIUM BETWEEN A TERNARY 
COMPOUND PHASE 


LIQUID AND A BINARY INTERMEDIATE 


An expression describing the equilibrium between a ternary liquid 
and a binary compound relates the heat of fusion and melting 
temperature of the binary compound, and the regular solution 
parameters, composition, and temperature of the ternary liquid 
phase. Following the steps outlined for the equivalent binary case 
given above, the following equation results. y 


&H c x ( 1 - T/T c ) 


= RxTx(N lc xln(N lc ) - N 2C xln (N 2C ) ) - RxTx (N, c xln (w, 
+ N lc xN 2C xL12 - N ic x(L12xN 2L x{1 - N 1L ) + ZUlxA/'-.x (1 

- L23xN 2L xN 1l ) - N ic x(L12xN 1l x (1 - N ) + L23xN-..x ( 1 

- L13xN 1l xN, l ) . 11 


+ N 2C xln(N 2L ) ) 
- N 1L ) 


(5) 


Equation (5) was solved as a function of temperature to yield the 
compositions of the Ga-As-Ge liquidus boundary which are in 
equilibrium with the compound GaAs . Figure (4a) shows the 

d ^ a ^f ain .^ nd the calculated liquidus boundary, which is 
the dominant liquid surface for this ternary system. 

This liquidus surface has two intersections with the previously 

boundar y' wh ich is in equilibrium with 
the Ge solid solubility phase, at a temperature of 873K. These two 
intersections define the compositional endpoints for which the 
liquid phase will be present. Figure 5a is a plot of the 1000K 
sectlon and shows the liquid and solid boundaries with 

COnn r tlng the liquid and solid phases at the two 
fan h!f tl0n endpoints mentioned above. Ternary compositions which 

tS^ 1 c-oi t iH een >, theSe ^ WO poin . ts Wl11 result in the presence of only 
two solid phases at equilibrium, the GaAs compound and the Ge 

r hSre i :i; ^ however ' a ternary eutectic involving GaAs, GeAs 
nf d th« e a Ge S '°i ld L solu . b ?- llt y Phase and will act to lower the limit 
wIth GaAs S .' r S solubility phase that is in equilibriim 

These calculations were also performed for the Ag-Ga-As ternary 
system The results of the previous section for the K 

2 b Xi? n '. Whlch determined the boundary between a ternary liquid 
phase that is in equilibrium with GaAs, is shown in Figure 4b. 

i re ® u J’ t ? this section were combined with the calculations 
reflected in Figure 4b and the resulting complete 800K Ag-Ga-As 
ternary diagram, showing liquidus and solidus boundarie^s with 
compositional tie-lines is presented in Figure 5b. 


TERNARY EUTECTICS OF THE Ag-Ga-As SYSTEM 

A eutectic point involving the ternary liquid phase, GaAs, GeAs 
and a primary Ge solid alloy phase was calculated from thi 
intersection of the three relevant ternary liquidus phase 
boundaries. These three boundaries are related to the equilibrium 
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between the ternary liquid phase and the GaAs compound, the GeAs 
compound, and the solid Ge alloy phase. The ternary liquid phase 
composition and temperature of the eutectic point and 
composition of the solid Ge alloy phase in equilibrium with the 
GaAs and GeAs at this temperature is listed in Table 4 . °| 

three solid phases in this eutectic reaction the Ge solid 
solubility phase contains 99.78 atomic percent Ge, 0.16 atomic 
percent As, and 0.06 atomic percent Ga. Ge compositions on the 
solidus boundary containing a higher percentage of As will be 
equilibrium with the GeAs compound. 

The composition and temperature of the ternary 

eauilibrium with the three compounds GaAs, GeAs, and GeAs 2 at the 
eutectic was determined from the intersection of the three liquidus 
surface boundaries described above but with the solid Ge alloy 
phase replaced by the compound GeAs 2 as the third liquidus boundary . 
These compositional values and eutectic temperature are also listed 
in Table 4 . 

Table 4 


Solid Phases 

Temp. (K) 

Liquidus 
Atomic % 
Ga As 

Balance 

GaAs, GeAs, Ge 

1002 

41.8 

1.7 

Gs 

GaAs, GeAs, 
GeAs 2 

1000 

60.9 

0.6 

G 0 


Calculated Eutectic Compositions for the As-Ga-Ge System 

Presented in Figure 6 is a representation of the complete As-Ga-Ge 
ternary phase diagram, including the three component binary 
diagrams which comprise the complete system. Visualization of the 
two eutectic points is enhanced by this representation. It is 
located by two distinct valleys which would form when the entire 
liquidus surface connecting the boundaries between the three binary 
diagrams is generated. 


SUMMARY 

Germanium and silver alloys compositions which are in equilibrium 
with gallium arsenide were determined from the appropriate ternary 
phase diagrams. The procedure for calculating these compositions 

consists of: 

1. Modeling the free energy functions of the binary and ternary 
phases by an application of the theory of a regular solution. 

2 Generating the representations for the binary and ternary 

liquid and solid boundaries from the calculated free energy 
functions . 
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r J gUlar solution parameters are calculated from the 
meltin g temperatures, and heats of fusion elemental values 
a nd the experimental values for the binary liquidus and 
sol idus boundaries . 

4. Approximating the ternary liquid and solid regular solution 
parameters as a sum of the component binary systems. 

5. The compositions of the ternary liquid and solid boundaries 
are calculated as a function of temperature. 

6. Calculating the compositions at the intersections of the 
ternary liquidus boundaries as a function of temperature. 
These compositions are the limiting values for which the 

P hase be P resen t . Between these limits only the 

solid phases will be present. ^ 

oht2?n t ^L tle ~ line f b6tWeen thS liquid and solid phases to 
obtain the compositions on the ternary solidus surface that 

are m equilibrium with the ternary liquid phase at ?he 
intersection points . 

8. The solid compositions resulting from Step 7 are the 

endpoints for which only two solid phases, GaAs and the 
terminal solid solubility phase will coexist. 

Ternary eutectic involving the ternary terminal solid solubility 

the com P° sltlona l range of solid alloys that will be in 
equilibrium with GaAs. a ternary eutectic point involvinq the 
three solid phases was obtained from the intersection of th? 

ern a ry liquidus surfaces whose compositions are in equilibrium 
with each of the three solid phases. equilibrium 

Experimental verification of the correctness of the ternary ohase 
equilibrium model is currently being undertaken. Differential 
scanning calorimetry is the primary technique being utilized to 
verify the equilibrium reactions predicted. Additional analysis to 
9 THq r f in h composition and compound formation is to include TEM and 
SIMS techniques but have yet to be initiated. Initial calorimetry 

react ions^wil 1 S ha ve>° indl< r ate that the non-equilibrium solid state 
reactions will have equal importance in the development of a 

suitable alloy system for high temperature metallizations to GaAs 
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FIGURE 2 

Temperature constitution phase diagram for the Ga-As system. Curve 1 is the 
calculated liquldus boundary and curve 2 the experimental boundary. 



FIGURE 3a FIGURE 3b 

The Ge and Ag rich corners of their respective ternary phase diagrams. Tie-lines 
shown connect the liquid and solid compositions which are in equilibrium. 
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baflsGe 673. K 
Dat,afuLe: ler^.plt 


GaRsAa 800. K 
DotoF'-le: ter^.plL 



FIGURE 4a FIGURE 4b 

Ga-As-Ge and Ga-As-Ag phase diagrams showing the ternary liquid composition 
which are in equilibrium with GaAs. 



FIGURE 5a 


FIGURE 5b 


Plots for the two liquid boundary intersections of the Ga-As-Ge and Ga-As-Ag ternary 
systems. Tie-lines show the liquid and solid compostions in equilibrium. 
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REPRESENTATION OF THE Ga-As-Ge SYSTEM 
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Ternary Ga-As-Ge phase diagram showing the 3 component binary diagrams. 



